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ABSTRACT 


The  electron  beam  sustainer  method  of  generating  a stable 
high  pressure  gas  discharge  is  investigated  theoretically  and  experi- 
mentally. The  most  important  discharge  instability  in  the  non- self- 
sustaining  discharge  is  the  streamer;  this  is  discussed  and  a model 
proposed  for  formation  and  propagation.  Integration  of  streamer 
velocity  over  the  discharge  gap  predicts  arcing  times.  Experimental 
measurements  of  streamer  behaviour  in  nitrogen  discharge  have  been 
made,  and  these  are  compared  with  predictions.  Comparisons  in  other 
gas  mixtures  and  various  electron-beam  sustainer  discharge  devices 
have  also  been  made. 


-111- 


d 


pr-f  - - 

i' 

fy 

■iA 

« 

, TABLE  OF  CONTENTS 

f 

- I 

6 Section 


Page 


I. 

II. 


III. 


Abstract 

List  of  Illustrations 

INTRODUCTION 

THEORETICAL  DISCUSSION 

A.  Streamers 

1.  Streamer  Model 

2.  Streamer  Current 

3.  Streamer  Tip  Radius 

4.  Streamer  Temperature 

5.  Streamer  Radius 

6.  Streamer  /Heating  at  Tip 

7.  Streamer  Velocity 

8.  Growth  of  Streamers 

9.  Streamer  Scaling 

10.  Summary 

EXPERIMENTAL  RESULTS 

A.  Preliminary  Cases 

1.  Rods  Parallel  to  the  Optic  Axis 

2.  Rods  Perpendicular  to  the  Optic  Axis 

B.  Framing  Camera  Experiments 

C.  Discussion  of  Results 

1.  Comparison  with  Theory 

2.  Predictions  of  Arcing  Limits  for  Discharges 

3.  Effect  of  Impurities 

D.  Initiation  of  Streamers 

E.  Conclusion 

Appendix  A - Gas  Specific  Heats 
REFERENCES 


iii 

vii 

1 

5 

5 

8 

15 

18 

19 

21 

21 

25 

30 

32 

34 


37 

37 

37 

38 
43 
52 
52 
57 
59 
62 
66 
67 
73 


j 

t 


I 

i 


-V- 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

1 Schematic  Drawing  of  Electron-Beam  Sustainer  Discharge,  3 

Showing  the  Electrodes  Used  for  Streamer  Investigations 

Z Pulsed  Discharge  Streamers  6 

3 Flowing  Discharge  Streamers  7 

4 Model  for  Streamer  Propagation  9 

5 Equilibrium  Electron  Density  for  N-  and  N,  + ^ Na  as  13 

Function  of  Temperature,  at  1 atm 

6 Mean  Thermal  Conductivity  for  Air  14 

7 Solution  to  Eq.  (5)  for  T as  Function  of  C 16 

8 Ratio  of  Cap  Sphere  to  Cylinder  Radius  for  Streamers,  as  20 

Function  on  Aspect  Ratio 

9 Streamer  Length  vs  Time  33 

10  Schematic  of  the  Optical  Arrangement  Used  for  Photographing  35 

Streamers 

11  Streamer  Growth:  Three  Cathode  Rods  Parallel  to  Optic  Axis  39 

12a  Streamer  Growth:  Single  Cathode  Rod  Perpendicular  to  40 

Optic  Axis 

12b  Experimental  and  Theoretical  Streamer  Length  42 

13  Experimental  Apparatus  Including  the  High-Speed  Framing  44 

Camera  Used  to  Photograph  Streamers 

14  Framing  Camera  Photo  Sequence  Showing  Propagation  of  45 

a Cathode  Streamer  and  Resultant  Arcing  When  the  Streamer 
Reached  the  Anode 

15  Growth  of  the  Streamer  Shown  in  Fig.  14  47 

16  Dependence  of  Streamer  Growth  Rate  Upon  Discharge  Power  49 

Density  and  Current  Density 


-vli- 


•'M, 


-I 


' f 

! i 


I i 

[ ri 


Figure 


Page 


17 

Dependence  of  Streamer  Growth  Rate  Upon  Discharge 
Power  Density  and  Electric  Field 

50 

18 

Streamer  Transit  Time  vs  Power  Loading  in  Nitrogen, 
for  Various  Electrode  Separations  (D) 

53 

19 

Streamer  Transit  Time  for  Cavities  with  Various 
Separations  Predicted  for  200°K  and  300”K  N_  Discharge 
at  1 atm 

55 

20a 

Loading  Energy  at  Streamer  Transit  in  Discharge  at 

1 atm 

58 

20b 

Loading  Energy  at  Streamer  Transit  in  Discharge  at 

1/ 10  atm 

60 

21 

Streamer  Critical  Radius  and  Latent  Period  for 

Spheroidal  Protrusion 

64 

A-1 

CVCP  Enthalpy  of  Nitrogen 

69 

A-2 

Constant  Pressure  Enthalpy  for  CO2 

70 

»V111- 


1 


u 

t : 

i i 

i i 
i 


I 


1 


IP 


'-I 

■ .'i 
•f' 


If 

I "■ 

i H 


1.  INTRODUCTION 

Stable  high-density  plasma  discharges  have  many  applications  in 
plasma  physics.  An  especially  significant  example  is  the  use  of  these 
discharges  in  high-power  gas  lasers.  The  work  described  here  concerns 
investigation  of  the  external  ionization  electron-beam  sustainer  method  of 
producing  a high-density  discharge,  together  with  the  adveuitages  and  insta- 
bilities peculiar  to  that  system. 

Classically,  stable  plasma  discharges  are  produced  by  working  in 
the  region  where  diffusion  of  the  charged  species  to  the  walls  stabilizes  the 
discharge.  The  electron  temperature  in  the  discharge  is  raised  (by  adjust- 
ing the  electric  field,  E/P)  until  the  ionization  rate  due  to  these  electrons 
equals  the  rate  at  which  electron-ion  pairs  recombine  at  the  walls.  How- 
ever, at  higher  pressures  and/or  in  larger  scale  devices,  the  diffusion  time 
generally  becomes  longer  thcun  other  instability  times  in  the  plasma,  and  the 
discharge  spokes  or  arcs.  Many  of  these  instabilities  are  associated  with 
the  ionization  process  and  its  coupling  with  the  electron  temperature.  AERL 
has  constructed  a facility  under  a company-sponsored  program  which  produces 
uniform  discharges  in  which  the  electron  density,  n^,  can  be  varied  independ- 
ently of  the  electron  temperature,  T^.  Separation  of  the  ionization  process 
from  the  electron  temperature  is  achieved  in  this  facility  by  having  a sepa- 
rate control  on  each,  resulting  in  a stable  discharge  over  a wide  range  of 
pressure  and  scale  size. 

We  achieve  the  separation  of  plasma  ionization  from  the  electron 
temperature  by  using  a beam  of  high-energy  electrons  (EB)  to  create  a 
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uniform  ionization  throughout  the  volume  of  the  plasma.  Using  the  high 
energy  electrons  to  provide  the  ionization,  the  electron  temperature  is 

varied  independently  by  adjusting  the  electric  field.  It  is  found  that  if  the  | 

electron  temperature  is  kept  so  low  that  the  ionization  due  to  the  discharge  j 

electrons  is  small  compared  to  the  ionization  due  to  the  high  energy  elec-  | 

1 

trons  the  discharge  remains  uniform  over  the  range  of  pressures  investi-  j 

i 

gated  (up  to  and  above  one  atmosphere).  As  the  electric  field  is  raised,  the 
ionization  by  the  discharge  electrons  becomes  substantial  and  the  discharge 
becomes  nonuniform  and  eventually  arcs.  Stable  discharge  behavior  is  ob- 
tained when  the  volumetric  electron  loss  rates  due  to  electron-ion  recom- 
bination and  electron-neutral  attachment  equal  the  volumetric  secondary 
electron  production  rate  by  the  EB.  The  recombination  coefficient  a and 
attachment  p are  important  parameters  determining  the  level  of  ionization 
obtainable  with  given  EB. 

Discharges  of  this  type  are  used  by  a number  of  ARPA-funded  lasers 
in  order  to  maximize  laser  output  from  high  pressure  gases.  In  addition  to 
10.  6jl  CO2  laser  discharges,  the  electron-beam  sustainer  principle  is  cur- 
rently being  used  for  CO  and  HF  lasers,  as  well  as  the  initiation  of  chemical 
fluorine-deuterium  lasers,  and  xenon  halide  visible  lasers.  In  view  of  the 
wide  range  of  application  of  this  discharge,  it  is  especially  valuable  to  in- 
vestigate its  stability  properties. 

The  discharge  apparatus  shown  in  Fig.  1 has  been  used  under  the 
present  contract  to  measure  electron  ion  recombination  and  electron  neutral 
attachment  rates  in  various  gases.  The  dissociative  attachment  instability 
discovered  under  this  contract  has  led  to  investigation  of  discharge  stability. 
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or  more  precisely,  modes  of  instability.  The  most  prominent  instability 
mode  in  the  externally  ionized  discharge  is  the  streamer,  and  some  progress 
has  been  made  in  understanding  this  phenomenon  in  the  present  case.  This 
will  be  discussed  below. 
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II.  THEORETICAL  DISCUSSION 
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A.  STREAMERS 

The  most  obvious  form  of  discharge  instability  occurring  in  non 
self-sustaining  discharges  is  the  streamer.  This  phenomenon  appears 
as  a visible  glowing  projection  from  one  discharge  electrode,  generally 
the  cathode  (although  anode  streamers  have  been  observed  under  certain 
circumstances).  The  projection  moves  across  the  discharge  at  a velocity 
near  that  of  sound,  and  when  it  reaches  the  opposite  side  an  arc  is  formed 
short-circuiting  the  discharge.  Figures  2 and  3 show  examples  of  streamers 
produced  in  an  electron-beam  sustainer  discharge.  Figure  2 is  a photo- 
graph of  a pulsed  discharge  in  the  gas  mixture  He;  N^:  CO^:  3:  2:  1 at 

12 

atmospheric  pressure.  With  electron  density  n - 3 x 10  cm'  and 

e 

1 2 

E/N  - 2 X 10  V cm  , the  cathode  is  at  the  lower  side  and  the  anode 
at  the  upper.  The  electron  beam  used  to  ionize  the  non- self- sustaining 
discharge  is  introduced  from  below,  and  the  electrons  are  of  such  a vol- 
tage that  they  do  not  densely  ionize  the  anode  region.  Since  current  must 
be  conserved,  the  electric  field  in  this  region  rises  sharply,  and  becomes 
high  enough  to  cause  the  generation  of  anode  streamers,  visible  at  the 
top  of  the  photograph.  Cathode  streamers  are  also  visible  at  the  bottom 
of  the  photograph.  Usually  the  electron  beam  energy  is  high  enough  so 
that  anode  streamers  hardly  exist  and  cathode  streamers,  growing  from 
discontinuities  associated  with  the  cathode  and  the  cathode  layer,  are  the 
most  important.  In  Fig.  3 the  streamers  of  a CW  flowing  gas  discharge 
are  shown.  Gas  flow  is  from  left  to  right  at  Mach  0.  4,  the  gas  used  is 
He:  N^:  CO^  16:  8:  1 at  380  torr  pressure,  and  the  cathode  is  at  the 

bottom,  from  which  the  electron  beam  enters  the  discharge.  Streamers 
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ANODE 


CATHODE  RODS 


Pulsed  Discharge  Streamers.  Electron  density  ~ 3 x 10  ( 

electric  field  5000  V/cm,  pressure  1 atm.  Gas  mixture 
He:  N^:  ^^2’  Electron  beam  enters  through  cathode 

bottom. 


are  visible  leaving  the  cathode:  they  then  are  blown  downstream  and  become 
complicated.  Anode  streamers  can  also  be  seen  leaving  the  upper  electrode, 
and  they  are  also  blown  downstream;  arcing  in  such  discharges  generally 
takes  place  several  cavity  widths  downstream,  presumably  when  the  dis- 
placed cathode  and  anode  streamers  finally  meet.  In  the  present  discus- 
sion we  will  consider  the  simpler  static  case. 

Attempts  to  explain  streamers  on  the  basis  of  electron  reproduc- 
tion rates  in  high- intensity  fields  have  generally  been  unsuccessful, 

(1,  V 

leading  to  predicted  velocities  ~ c/lOO,  where  c is  the  velocity  of  light. 
Although  these  velocities  are  correct  for  the  initiation  of  spark  breakdown 
electron  avalanche  in  gases  subjected  to  high  electric  fields,  this  type  of 
avalanche  cannot  occur  in  the  low  fields  of  the  nonselfsustaining  discharge. 
The  process  forming  streamers  in  such  a discharge  must  be  a gasdynamic 
heating  effect  in  which  gas  near  a discontinuity  is  heated  to  a sufficiently 
high  temperature  so  that  it  becomes  conductive;  the  discontinuity  is  then 
effectively  propagated  out  through  the  discharge  in  the  form  of  a hot,  highly 
conducting  column  of  gas.  We  will  examine  the  properties  of  such  a 
streamer,  simplifying  the  discussion  as  much  as  possible  in  order  to 
obtain  analytical  and  avoid  numerical  solutions. 

1.  Streamer  Model 

We  will  regard  the  streamer  as  a cylindrical  column  of  gas  at 
temperature  T,  with  radius  b and  length  c.  We  approximate  the  head  of 
the  streamer  as  a sphere  of  highly  conducting  gas  of  radius  a,  in  which 
joule  heating  of  the  gas  occurs.  (See  Fig.  4)  The  pressure  inside  the 
streamer  is  assumed  identical  to  that  in  the  main  discharge  (if  it  were 
higher  the  streamer  would  expand  in  all  directions  near  sonic  velocity 
into  the  discharge  gas). 


The  volumetric  joule  heating  in  the  cylinder  will  then  be  = — 

(Trb'^j'^CT 

where  o is  the  gas  conductivity  and  I is  the  current  carried  in  the  cylinder. 


^4  o. 


For  temperatures  below  - 3 x 10  K,  in  atmospheric  pressure  gas  conductive 
heat  loss  will  greatly  exceed  radiative  loss.  The  thermal  conduction 
equation  for  a cylinder  may  be  written  as 


ZTrr/c^-^^  = -Trr^q,  where  q is  the  volumetric  heating  rate,  so  that 


. 2 dT 

~ ■ R dR  ■ 


(3) 


This  may  be  approximated,  following  Raizer,  as 

T 


f KdT. 
R 


(1) 


0 


where  A is  a numerical  factor  and  the  integration  is  a convenient  way  of 


using  a mean  value  of  the  conductivity,  since  this  quantity  varies  greatly 


.4  o. 


over  the  temperature  range  considered  (T  <10  K). 

T 


Writing  K ^ ^ J K dT  , and  noting  that  for  T <3x10^  radiation 
losses  are  negligible  compared  to  conduction  losses,  we  equate  the  heating 
and  cooling  rates  and  obtain,  placing  R = b. 


,2 

7rb  a 


= ttAIcT 


(2) 


Adopting  Raizer 's  value  for  high  temperature  air,  A ft  15,  we  write  the 
conductivity  as  o = n^  e|i^/N,  where  n^  is  electron  density,  e is  electron 
charge,  is  electron  mobility  and  N is  gas  density.  We  will  determine 


n from  the  solution  for  the  Saha  equation  in  the  gas  at  the  temperature 
e 
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T considered,  and  the  density  will  be  inversely  proportional  to  temperature 
under  the  assumption  of  constant  pressure,  so  that  the  conductivity  scales 


and  we  obtain 


• w 


where  a , T and  n are  the  values  for  the  electrical  conductivity, 
temperature  and  electron  density  in  the  undisturbed  medium,  n^(T)  is 
obtained  from  the  appropriate  Saha  equation,  and  n^(T^)  « Here 

n^^  corresponds  to  the  density  of  electrons  produced  by  external  ioniza- 
tion. The  dissociation  of  molecular  gas,  which  will  occur  at  the  tempera- 
tures attained  in  the  streamer,  will  increase  the  electron  mobility  and  the 
conductivity.  In  Eq.  (3)  we  have  ignored  this  effect  since  the  change  in 
n^(T)  with  temperature  is  much  more  rapid  and  has  much  the  greater 
effect  on  the  conductivity.  Equation  (3)  now  expresses  the  relationship 
between  the  current  in  a cylindrical  gas  column  of  radius  b and  its 
temperature,  in  the  regime  where  the  dominant  cooling  process  is  con- 
ductive. In  order  to  use  this  we  need  the  function  n^(T)  from  the  Saha 
equation  and  the  current  I,  The  first  of  these  is  given  by  the  expression 


_^o^o^o  t r«l,  . 4 .1 
% T 2 L’*  - ‘J  • 


where 


f -* . 


^ u T L p l~r5  I 
o o o*^  \ h / 


27rm^ky/2  -i/kT 


(4) 


and  L , T and  P are  Loschmidt’s  number,  initial  gas  temperature 
o o o 

(°K)  and  pressure  in  atmospheres,  u^  and  u^  are  the  statistical  weights 

of  the  ion  and  atom,  p is  the  degree  of  dissociation  for  a molecular  gas 

(P  = 2 for  nitrogen  in  the  regime  of  interest),  m^,  k and  h are  the 

electron  mass,  Boltzman  constant  and  Planck  constant  respectively, 

(4) 

and  I is  the  ionization  potential.  Values  of  n^  derived  for  pure  N^, 

and  for  contaminated  with  a molar  fraction  4*  of  Na,  are  shown  in 

Fig.  5 at  atmospheric  pressure.  ^ 

The  mean  thermal  conductivity,  J KdT,  is  derived  from  the 

o 

published  conductivity  for  air^^^  at  high  temperatures,  which  is  close  to 

that  for  nitrogen,  and  is  shown  in  Fig.  6.  Over  the  range  of  interest 

this  quantity  can  be  fitted  fairly  well  by  the  expression  k = ‘ ^^^w/cm/°K, 

where  k =3.  145  x lO"^^  and  T is  in  degrees  Kelvin, 
o 

Combining  Eqs.  (3)  and  (4)  with  the  above  expression  for  conductivity 
yields  an  implicit  equation  for  the  gas  temperature  at  which  the  streamer 
can  exist; 


f«( 

I 

i tl 

I 


The  solution  of  Eq.  (5)  for  10^^  S C < 10^^  is  shown  in  Fig.  7 for  nitro- 
gen. It  is  evident  that  the  streamer  temperature  is  a slowly  varying  func- 
tion of  C,  and  5000  < T < 10,  000  over  this  range  of  C.  We  note  in  passing 
that  if  the  Saha  equation  is  solved  for  gas  containing  an  impurity  of  low 
ionization  potential,  such  as  sodium,  the  stable  streamer  temperature  is 
considerably  reduced,  and  solutions  to  Eq.  (5)  are  shown  for  nitrogen 
with  added  sodium  contaminant  of  molar  fraction  4^,  We  might  therefore 
anticipate  faster  formation  and  propagation  of  streamers  in  contaminated 
gas. 

2.  Streamer  Current 

The  next  problem  is  to  estimate  the  current  I in  the  streamer. 

Most  of  the  current  will  be  collected  in  the  high-field  spherical  region 
at  the  tip  of  the  streamer.  In  order  to  calculate  this  we  need  to  know  the 
potential  distribution  near  the  tip.  We  will  therefore  make  the  following 
assumptions: 

(1)  The  field  can  be  approximated  by  that  of  a sphere  of  radius 
a between  two  plates:  denote  the  sphere  potential  V^,  the 
top  plate  V^,  bottom  plate  zero  and  plate  separation  D,  dis- 
tance from  sphere  to  bottom  plate  c. 

(2)  All  current  enters  the  streamer  through  the  upper  half  of 
the  sphere. 

(3)  The  sphere  is  attached  by  a thin  conducting  cylinder  (the 
streamer)  to  the  lower  plate. 

(4)  The  field  is  E = V_/D  at  infinity. 

o t- 

The  continuity  equation  then  gives 

V • j = q = 0 
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where  j is  current  density  and  q is  charge  density.  This  may  be  written 


V • (aE)  = 0 


Vct  • VV  + (jV‘"V  = 0 (6 

Since  we  are  approximating  the  streamer  tip  as  a conducting  sphere, 
the  conductivity  will  be  constant  outside  the  sphere  and  equal  to  the 
undisturbed  value  o^,  and  Eq.  (6)  becomes  Laplace's  equation 

V = 0 (7 

The  potential  field  for  a charged  sphere  in  an  infinite  field  must 


then  be  of  the  form; 


^1  / 

■ 1^2  ^ ^ 


using  polar  coordinates  (r,  0)  centered  on  the  sphere,  higher  terms  being 
eliminated  by  the  boundary  condition  (4).  This  expression  is  slightly 
erroneous  near  the  electrode  plates,  where  image  effects  increase  the 
field,  but  the  error  is  xinimportant  for  most  of  the  streamer  motion 
between  the  plates,  and  will  be  neglected.  Noting  that  a^  = “^2/^ 
the  sphere,  and  = a^  + ^^1^,  we  eliminate  the  constants  and  obtain; 


^"V-  a7c  ’^l  cd-aVc^  + (a/r){Vj.E^c(l-aVc^)}] 


Differentiating  Eq.  (8)  and  placing  r = a we  obtain  the  radial  field  at  the 
surface  of  the  sphere,  and  since  under  out  assumptions  the  sphere  will 
be  at  virtually  the  same  potential  as  the  lower  plate,  we  take  Vj  = 0 and 


obtain: 


Er=a'E„{[l+f  +f]  + 3co,4 


The  current  entering  the  sphere  will  then  be 


n 

= I E 

J o r-a 


Ztt  a sin  Q d0 


I = 4 TT  a a E 
o I 


o B ^ f] 


where  the  second  term  in  Eq.  (9)  has  been  integrated  over  the  upper  hemis- 
phere, and  since  — is  generally  greater  than  10,  the  first  term  in  Eq.  (9^ 
contributes  most  to  the  current;  only  a small  error  is  introduced  by 
ignoring 

3.  Streamer  Tip  Radius 

We  wish  to  estimate  a proper  radius  for  the  sphere  at  the  tip  of 
the  streamer.  A reasonable  assumption  is  that  this  radius  corresponds 
to  the  value  at  which  the  electric  field  is  high  enough  for  the  secondary 
electrons  of  the  discharge  to  produce  further  ionization.  While  the 
increased  conductivity  of  the  sphere  caused  by  this  further  ionization 
will  reduce  the  field,  the  value  obtained  from  the  solution  to  Laplace's 
equation  will  be  approximately  correct.  We  therefore  write  the  field 
intensification  a at  distance  r from  the  center  of  the  hemispherical 
cap  on  the  streamer,  radius  b and  height  c,  by  differentiating  Eq.  (8)  as: 
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(10a) 


a = 


cos  6, 


where  0 ~ 0,  since  the  solution  ignores  the  field  from  the  streamer 
cylinder  and  considers  only  that  near  the  tip  due  to  the  hemisphere. 

r 

Taking  0 = 0 and  solving  this  equation  for  we  obtain  the  values  shown 
I.  Pig.  8 for  various  values  of  the  intensification.  Typically  ionization 
IS  important  for  - 3E^,  and  a representative  value  is  ^ = 3.  Identify- 
ing the  radius  with  the  cap  sphere  radius  a,  we  may  express  this  quantity 
from  Eq.  (10a)  in  good  approximation  as: 


a^ 

'b 


V 


1 + c/b~ 
a - 1 ’ 


where  b/c  is  neglected, 


and  since  c/b  » 1,  we  obtain 


^ te^c/b  (q  - 1) 


(11) 


4.  Streamer  Temperature 

We  can  now  combine  Eqs.  (5)  and  (10)  and  use  the  above  value  for 
a/b  to  estimate  the  temperature  of  a thermally  stable  streamer  column. 
The  expression  for  the  constant  becomes: 


C = 


U 

15 


T (j  E ^ 
o o o 


K 


o 


(12) 


We  note  that  C depends  only  slightly  on  the  streamer  radius.  For  a 
typical  discharge  in  nitrogen,  we  have  c - 1 cm;  taking  b = 1 mm, 

T^  = 300°K,  = Z.  6 X lO”^  mho  cm  \ E^  = 5000  v/cm,  and  recalling 
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that  = 3.  145  X 10  we  get  C = 6.  4 x 10^^.  From  Fig,  7,  this 
corresponds  to  a tennperature  T = 6200°K  in  pure  nitrogen. 

Thus  we  have  obtained  the  temperature  at  which  the  streamer 
column  will  be  stable:  we  note  that  it  is  positively  stable  since  a local 
increase  in  temperature  will  tend  to  increase  electrical  conductivity 
and  reduce  joule  heating,  as  well  as  increase  the  cooling,  so  that  the 
temperature  will  be  reduced;  similarly  a local  decrease  in  temperature 
has  the  opposite  effect. 

5.  Streamer  Radius 

The  radius  b of  the  streamer  is  not  determined  by  stability 
criteria,  which  merely  fix  the  temperature  with  respect  to  current 
passing  through  the  streamer.  We  will  set  a limit  on  the  radius  by 
st^posing  that  the  streamer  expands  radially  at  a velocity  corresponding 
to  the  rate  of  heat  diffusion  into  the  ambient  gas.  The  visible  streamer 
will  be  the  gas  at  temperature  near  6000°K,  which  will  also  be  conducting 
most  of  the  current.  For  discharge  duration  T,  the  streamer  base 
radius  at  the  end  of  the  discharge  will  then  be  b ft  ^/c  ^/pCp.  We  note 
that  b will  increase  as  discharge  pressure  decreases.  Taking  the  values 
corresponding  to  10^  °K  for  we  obtain  b ft  ^4.5  t/N,  where  N is 
gas  density  (amagat),  and  typically  in  an  atmospheric  discharge  after 
100  jjisec  we  obtain  b ss  0.  2 mm.  We  expect  streamers  to  be  narrower 
in  the  shorter,  more  highly  pumped  discharges. 

6.  Streamer  Heating  at  Tip 

We  will  now  estimate  the  rate  of  heating  of  the  gas  in  the  region 
of  the  streamer  tip.  As  this  gas  is  heated  it  will  expand  isotropically 
thus  increasing  the  electrical  conductivity,  and  at  high  enough  temperatures 


-?I- 


sufficient  thermal  ionization  will  be  present  to  further  increase  the  con- 
ductivity. 


Since  the  electrical  conductivity  must  rise  very  rapidly  as  r de- 
creases in  the  sphere,  an  appropriate  model  for  the  conductivity  might  be 
ka/r 

a = ' e ' where  k is  a constant  and  the  conductivity  at  r = a is  the 

- k 

same  as  that  in  the  undisturbed  medium,  so  o ' = a e . The  net 

’ o o 

• z 

heating  rate  q within  the  sphere  is  then  obtained  by  integrating  j /a 
over  the  cap  sphere,  down  to  a radius  corresponding  to  the  cylinder 
radius.  We  obtain: 


« - j 2,2  -r-r-  '>■■■ 

^ (47rr  ) o 


This  is  integrated  by  substituting  p = — , and  we  obtain 


-k(a/b  - 


• watts 


Noting  that  the  ratio  of  conductivities  at  6200°K  and  300°K  is  approxi- 
3 

mately  10  , allowing  for  the  external  ionization,  and  taking  the  value 
a/b  = 3,  we  obtain  k = 3.  45.  Thus  the  second  term  in  the  expression 
for  q can  be  ignored,  and  the  mean  volumetric  heating  rate  over  the 
sphere  is: 


3r  / 3 

J—  • W/cm 


16  TT  0 a k 
o 


This  result  does  not  depend  strongly  on  the  model  chosen  for  the 
conductivity,  for  if  instead  we  write  o = 0^(r/a)  where  m is  a constant. 


in  order  to  represent  the  rapidly  increasing  conductivity  as  the  radius 

2 

diminishes,  we  can  integrate  j /cj  over  the  sphere  and  average  it  as 
before,  and  we  obtain: 


2 4 

1 6 TT  CTq  a (m-  1) 


in  place  of  Eq.  (13).  For  the  same  range  of  temperatures,  conductivities 
and  radii  we  find  m = 6.  29;  the  ratio  of  the  two  volumetric  heating  rates 
is  then  3/2  approximately,  so  that  the  determination  of  is  close  to 
model- independent.  In  the  subsequent  work  we  will  use  that  given  by 
Eq.  (13). 

We  have  ignored  radial  temperature  differences  due  to  variation 
of  heating  rate  radially  within  the  sphere;  in  actual  fact  only  the  central 
region  will  reach  the  cylinder  temperature  derived  in  Eq.  (5),  and  we 
are  over-estimating  the  energy  required  to  do  so.  An  approximate  esti- 
mate of  the  magnitude  of  the  correction  for  the  fact  that  the  temperature 
varies  inversely  with  radius  in  the  sphere  may  be  made  as  follows. 


k(a/r-  1 ) 

We  have  written  the  electrical  conductivity  as  a = e ' 'in  Eq.  (12), 

and  shown  that  the  volumetric  heating  rate  is  relatively  model- independent. 

It  is  a good  approximation  for  the  temperature  dependence  of  conductivity 

T/T^  o 

in  the  region  T > 6000°  to  write  a = where  898  K 

and  a base  conductivity  is  provided  by  external  ionization.  The  self- 
consistent  variation  of  temperature  with  radius  is  then 


T = k T,  (a/r-1). 
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The  volumetric  heating^needed  to  bring  the  center  of  the  sphere  to  the 
critical  temperature  T at  which  it  becomes  sufficiently  conductive  to 
support  the  streamer  cylinder  is  then; 


T 


<T1>  = ^ j 47r  r^  (T/Tj)^dr 

b 

where  V is  the  sphere  volume  and  Y = 0.  63  for  a nitrogen  discharge 
(see  Appendix  A).  Making  the  substitution  and  comparing  this  value  with 
the  volumetric  heating  t]  needed  to  bring  the  entire  sphere  to  the  tempera- 
ture T,  we  obtain  the  correction  factor 

= I(b/a), 

where 

1 

I(b/a)  = j x^(l/x-l)^  dx. 
b/a 

Numerical  integration  of  I taking  Y = 0.  63  is  given  in  Table  I, 
showing  that  I is  a slowly  varying  function  of  b/a,  and  for  b/a  = l/3 
we  obtain  1 = 0,  158.  For  k = 3.45,  T2  = 898°K,  T = 6200°K,  we  have 
^h^/h  = 0.30.  We  conclude  that  ignoring  the  temperature  distribution 
in  the  sphere  will  result  in  calculated  heating  that  is  too  large  and  should 
be  reduced  by  a correction  factor  X = 0,  3, 


1 


TABLE  I 

NUMERICAL  INTEGRATION 


b/a 

I 

X = <Ti  >/r\ 

0,  180 

0.  175 

0.  3 

0.  163 

0,  316 

0.  333 

0.  158 

0.  306 

0.4 

0.  146 

0,  283 

0.  5 

0.  124 

0.  240 

7.  Streamer  Velocity 

At  this  point  we  can  identify  a characteristic  velocity  for  propa- 
gation of  the  streamer.  It  is  the  radius  of  the  sphere  divided  by  the  time 
taken  to  reach  the  stable  streamer  temperature  obtained  above.  Using 
the  Air  Force  values  for  the  enthalpy  of  N_  (see  Appendix  A),  we  can 
write  the  energy  required  to  bring  1 cc  of  gas  at  1 atm  from  300°K  to 
temperature  T as: 

e =6^(T/t^)°-^^ 

where  T^  = 3500°K,  = 1 joule/cm  . In  the  present  model,  the 

mass  of  gas  heated  is  continuously  diminished  as  it  expands  outward 
past  r = a and  leaves  the  region  of  heating;  we  can  ignore  the  conductive 
heat  loss  compared  to  this  large  convective  heat  loss,  which  is  auto- 
matically included  in  the  expression  for  €,  and  the  total  volumetric 
heat  input  under  these  conditions  will  be: 

Q = eNX  (14) 

where  N is  gas  density  normalized  to  1 atmosphere,  taking  initial 
temperature  as  300°K. 
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J 


I 


I 
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Identifying  this  energy  with  the  heat  input  gives  the  heating  time  T: 


^ = Q/qs- 


Writing  the  streamer  velocity  as  v = a/r  and  collecting  terms  we  obtain 
the  velocity  as; 


V 


3 a o E ^ 
o o 

X k e N 


(15) 


For  a typical  discharge  in  nitrogen,  taking  the  values  a = 0.  4 mm, 
Oq  = 2.  6 X 10  ^ mho  cm  = 3500  V/cm,  c = 0.  5 cm,  X = 0.  3, 

= 1.  12  X 10"^  g/cm^,  = 3500°K,  T = 6200°K,  k = 3.45  we  get 

4 / 

V = 1 . 5 X 1 0 cm/ sec 


and  we  note  that  the  sound  speed  in  the  undisturbed  gas  is 

V = 3.  53  X 10^  cm/sec,  so  that  in  this  case  the  streamer  travels  at  less 
c 

than  Mach  0.  5,  implying  that  the  assumption  of  uniform  pressure  is  still 
consistent  with  the  predicted  velocity. 

For  velocities  below  Mach  1 the  assumption  that  the  gas  expands 
rapidly  enough  so  that  the  expression  for  the  energy  input  per  unit 


volume  required  to  reach  temperature  T,  € 


1 


where  gas 


o 

density  p is  a function  of  T,  will  be  approximately  valid.  In  this  region 
the  constant  volume  enthalpy  discussed  in  Appendix  A should  be  used  in 
calculating  the  velocity  from  Eq.  (15). 

At  supersonic  speed,  however,  the  gas  cannot  expand  rapidly  enough 
during  the  heating  process,  and  the  gas  density  will  remain  close  to  its 
initial  value  p^  during  heating.  The  energy  input  per  unit  volume  is  then 
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^3  =Po 


dT,  and  the  streamer  velocity  must  be  calculated  based 


on  the  internal  energy.  Since  the  density  is  still  at  the  initial  value, 

and  the  discharge  power  loading  into  the  streamer  must  be  cor- 
respondingly greater  in  order  to  make  the  streamer  propagate  supersonically. 

In  the  intermediate  region  the  energy  input  per  unit  volume  re- 
quired to  reach  temperature  T has  an  intermediate  value,  or  g <€o- 

While  we  do  not  know  exactly  how  this  energy  is  distributed  between  gas 
internal  energy  and  gas  kinetic  energy,  we  do  know  that  the  resultant 
streamer  velocity  must  be  Mach  1.  Thus  we  can  distinguish  subsonic 
sonic  and  supersonic  regimes  of  streamer  propagation.  Since  the  power 
loading  into  the  streamer  increases  as  it  moves  across  the  discharge, 
we  expect  these  regimes  to  succeed  each  other  in  sequence. 

The  subsonic  regime  is  considered  in  the  discussion  above  leading 
up  to  Eq.  (15).  This  equation  will  be  taken  as  valid  with  g = until  the 
velocity  equals  the  soundspeed  in  the  undisturbed  medium,  a point  we 
identify  as  c = c^.  After  this  point  the  velocity  will  be  assumed  to 
remain  at  the  soundspeed  until  the  streamer  has  reached  a length  c = c^, 
at  which  the  velocity  begins  to  become  supersonic.  As  the  streamer  fur- 
ther increases  in  length  the  velocity  will  continue  to  increase.  Since  the 
only  physics  altered  is  the  value  of  £,  Eq.  (15)  can  still  be  used  to  derive 
the  supersonic  velocity,  placing  € = There  is  a small  change  to  be 

made  in  the  value  of  in  the  case  of  nitrogen  at  1 amagat  the  internal  energy 

T/T^  _3 

may  be  suitably  approximated  as  6^  = e , where  6^  = 0.  596  j cm 
and  T^  = 2010°K.  We  may  then  derive  \ as  before  and  obtain 
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I 


X 


3 


3e 


.(T  + kT2)/T4 


I 


2 

X 


k T /(T  X) 
e dx 


(16a) 


where  k = 3.45,  T = 6200°K  and  = 898°K,  and  we  define  the 

supersonic  value  of  X.  Taking  b fe  a/3,  as  we  did  in  the  previous  case, 
and  evaluating  Eq.  (I6a)  we  obtain  X^  * 0-  105,  which  will  be  used  below. 

We  may  now  use  Eq.  (15)  to  determine  the  boundaries  c^  and  c^  of 
the  sonic  regime.  In  order  to  do  this  we  must  identify  the  value  of  the 
core  radius  b in  the  head  region  of  the  streamer.  We  will  assume  that 
this  is  the  propagation  distance  of  the  conduction  heat  wave  during  the  time 
the  discharge  head  moves  through  its  radius,  keeping  in  mind  the  fact 
that  the  gas  in  the  central  region  of  the  cap  sphere  will  be  at  the  stable 
streamer  temperature,  6200°K  for  N^.  Thus 


b = ^ f where  V is  the  streamer  velocity. 


K a 


pC 


(16b) 


Combining  this  with  Eq.  (11)  we  obtain  the  expression 


a = f c 


2/3 


(16c) 


rx/pC 

where  f = , and  substituting  in  Eq.  (15)  the  velocity  may  be 

L(a-l)^  J 


written 


dc  _ 2/3  \ , ^i 

^-A.  /Ci  [V4+—  J 


(17a) 
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Vi 


where  = 3ct  ^ €•)•  and  i = 1 or  3 for  the  subsonic  and  super- 

sonic limits  of  the  sonic  regime,  respectively.  For  a nitrogen  discharge 
we  take  f = 0.  055  [^since  V is  present  to  the  l/3  power  it  is  a slowly 


varying  factor  and  we  can  take  V = soundspeedj,  and  identifying  with 


the  soundspeed  we  can  solve  Eq.  (17a)  for  the  limits  c^  and  c^.  These 


are  given  in  Table  II  for  a nitrogen  discharge  at  atmospheric  pressure 


and  300°K.  For  typical  loadings  near  3 kW/cm^,  the  streamer  will 


move  at  sonic  velocity  between  0.  3 and  1. 4 cm,  thereafter  becoming 
supersonic.  Streamer  behaviour  as  a function  of  time  will  be  discussed 
in  the  next  section. 


TABLE  II 

STREAMER  SONIC  LIMITS 


cr  E ^ 
o o 

W cm"^ 

‘^l 

cm 

^3 

cm 

10^ 

0.  84 

3.  43 

3x10^ 

0.  34 

1.  39 

lo" 

0.  134 

0.  57 

3 X lo'^ 

0.  06 

0.  24 

10^ 

0.  02 

0.  09 

Finally  we  notice  that  the  initial  value  of  b at  the  electrode,  b^, 


dc 


may  be  estimated  by  taking  b^  = c and  V = combining  Eqs.  (16b),  (16c) 
and  (I7a)  yields  the  transcendental  equation: 


1 


r 


>. 

c 

p*- 

\ , 


(17b) 


For  an  atmospheric  pressure  discharge  we  find  the  values  for 
in  Table  III.  We  take  as  the  value  corresponding  to  a typical  discharge 
b =0.016  cm. 

° TABLE  III 

STREAMER  RADIUS 


r 2 

(j  E 

0 0 

W cm  ^ 

b 

0 

cm 

3 

^ -2 

10 

2.  9 X 10 

-2 

3x10 

1.8x10 

4 

-2 

10^ 

1.0x10 

5 

-3 

10 

3.  4 X 10 

8.  Growth  of  Streamers 

The  length  of  a streamer  for  c < or  for  c > C2  will  be  given 
by  integration  of  Eq.  (17a)  between  appropriate  limits,  writing  c^  = c. 

In  the  intermediate  region  the  streamer  will  move  at  the  soundspeed.  We 
consider  the  three  regimes  separately. 

a)  Subsonic  Regime 

The  lower  integration  limit  is  the  minimum  radius  at  the  electrode, 
b^,  and  the  upper  is  distance  from  the  electrode,  with  the  condition 
Dj  ^ c 2 • have 


i 

I 


I 


i 

i 
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<■( 


/ 


c-2/3  r7  .c'/M-" 


J dc=A,Tj, 


where  is  elapsed  time.  This  becomes: 


(18a) 


r,  b'/h' 


1 - (b^/D,)‘/^ 

1 + 7f /(4D 


(18b) 


The  transition  to  the  sonic  regime  occurs  at  the  limit  = c^ 
b)  Sonic  Regime 

In  the  intermediate  regime  the  velocity  is  sonic,  and 

r2  = (D2-<=i)/\’ 

where  is  distance  from  electrode  and  is  soundspeed,  and 

Cl  <D2  <03. 


(I8cl 


The  integration  limits 

are 

now  c 

/■°3  ^-2/3 

7 + 

C-V3] 

J f 

4 ^ 

f 

^7. 

dc  = A3  T3 


(18d) 


Integration  yields: 


T — 

f 


1- 

1 + 


-(C3/D3) 


1/3 


7f/(4D3^/^ 


(18e) 


I 


with  the  condition  D3  > C3. 
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Eqviations  (18b),  (18c)  and  (18e)  now  give  us  the  length  of  time  taken 
for  the  streamer  to  reach  any  given  length.  Formally  speaking  it  can  be 
see  from  (18e)  that  the  streamer  will  cross  any  discharge  in  a finite  time; 
letting  Dj -*  oo  yields  Tj -►  3/[A^(7/4  The  history  of 

streamer  length  in  an  atmospheric  pressure  nitrogen  discharge  is  shown  as 
function  of  elapsed  time  in  Fig.  9 for  discharges  at  various  power  loadings. 

It  is  apparent  that  the  sonic  regime  occupies  a relatively  constant  fraction 
of  the  total  duration.  We  would  expect  streamers  to  traverse  relatively 
large  discharges  30  cm  across  almost  as  rapidly  as  those  only  5 or  10  cm 
across.  Since  the  discharge  arc  initiates  according  to  the  present  hypothesis 
when  the  streamer  has  traversed  the  discharge  cavity,  we  expect  the  ad- 
vantage in  arcing  limit  gained  by  large  discharge  electrode  spacing  to  be 
very  small.  Predictions  of  arcing  limits  and  comparison  with  experiment 
will  be  given  below, 

9.  Streamer  Scaling 

It  is  interesting  to  note  that  in  the  subsonic  regime  the  transit 
time  given  by  Eq,  (18b)  is  essentially  the  constant-pressure  energy  re- 
quired to  heat  the  gas  up  to  temperature  T given  by  Eq.  (5)  divided  by 
the  rate  of  heat  input.  If  E/N  is  kept  constant  and  the  gas  density  N is 
varied,  for  most  applications  of  the  non- self- sustaining  discharge  it 
is  kinetically  favorable  to  maintain  n^^/N  constant,  which  results  in 
unchanged  pumping  rate  per  molecule.  The  electrical  conductivity  will 
then  be  independent  of  gas  density,  and  T scales  as 

t = K S ^ (19a) 

a (E  /N)  N 
o'  o'  ' 
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so  that  arcing  time  varies  inversely  with  gas  density,  which  is  consistent 
with  the  observation  thct  EB  sustainer  discharges  in  gases  at  low  pressures 
arc  less  readily  than  at  high  pressure  for  the  same  E/N.  The  same  scaling 
applies  for  the  time  spent  in  the  supersonic  regime. 

In  practical  discharges  the  streamers  are  likely  to  reach  the 
sonic  or  supersonic  regimes  before  completing  their  journey.  We  note 
that  since  the  soundspeed  is  independent  of  gas  density  low  pressure  dis- 
charges of  the  CW  type  are  less  likely  to  attain  the  supersonic  regime, 
and  therefore  will  take  longer  to  arc,  than  high  pressure  discharges. 

The  values  of  c^  and  have  been  computed  from  Eq.  (17a)  for  the  gas 
mixture  He:  N2:  CO2  3:  2:  1 at  various  pressures  and  power  loadings, 
and  are  shown  in  Fig.  10,  Equation  (17a)  reduces  fairly  accurately  to  the 


r . 3/4 


, since  f = N 


expressio 

duration  of  transit  of  the  sonic  regime  varies  as 


-1/3 


, and  thus  the 


T 


2 


(19b) 


Consequently  the  lower  the  pressure,  the  longer  the  relative  duration  in  the 

2 

sonic  regime.  The  variation  of  transit  time  T with  power  loading  o^  E^ 
in  the  sonic  case  is  not  quite  as  fast  as  in  the  subsonic  and  supersonic 
cases,  both  of  which  follow  Eq.  (19a),  thus  leading  to  slight  variations 
from  the  straightforward  concept  of  a limiting  permissible  discharge 
loading  energy.  Predictions  will  be  made  for  various  discharges  below. 

10,  Summa  ry 


The  streamer  model  developed  above  derives  a streamer  velocity 
proportional  to  the  square  of  the  discharge  electric  field  in  the  subsonic 


Fig.  10  Schematic  of  the  Optical  Arrangement  Used  for  Photographing 
Streamers 


regime.  After  a sonic  transition,  at  high  fields  the  predicted  velocity- 
will  be  supersonic  and  the  assumption  of  uniform  pressure  will  be  invalid, 
resulting  in  corrections  to  the  thermal  equations.  The  electrical  con- 
ductivity is  an  extremely  steep  function  of  temperature,  and  as  soon  as 
the  central  region  of  the  sphere  reaches  sufficient  temperature  the  ioni- 
zation region  moves  on,  ingesting  fresh  gas.  Only  the  core  reaches  the 


high  temperature  of  the  luminous  streamer. 

The  streamer  cylinder  will  expand  radially  at  a velocity  corre- 
sponding to  the  thermal  diffusion  rate,  so  that  it  will  acquire  the  appearance 
of  a cone  the  tangent  of  whose  half  angle  will  equal  the  ratio  of  the  axial 


and  radial  velocities.  However,  taking  the  radius  b = 0.  05  cm,  the 


radial  velocity  is  approximately  v fc 
is  much  less  than  the  axial  velocity. 


K / 

— >,■  T ft:  90  cm/sec  for  N_,  which 

pCpb  ' L 

and  the  above  derivation  of  thermal 


equilibrium  for  a cylinder  will  be  satisfactory.  Approximately  speaking, 
arcing  times  scale  inversely  with  power  loading  and  gas  pressure. 


III.  EXPEEIMENTAL  EESULTS 


A.  PEELIMINAEY  CASES 

An  electron-beam  sustainer  discharge  device  was  used  to  produce 

streamers  in  a nitrogen  discharge  at  atmospheric  pressure.  The  device, 

f7  8) 

which  has  been  described  previously,  ' ’ ^ is  shown  schematically  in  Fig.  1. 
A 90  kV  electron  beam  enters  the  chamber  through  a 1/2  mil  aluminum  foil 
and  an  electrode  plate  perforated  with  1/8  inch  holes  in  a 5 x 10  cm  rectan- 
gular array;  it  ionizes  the  gas  present  and  a discharge  then  takes  place  in 
the  secondary  electrons  between  the  perforated  plate  and  the  second  elec- 
trode; discharge  spacing  is  2.  54  cm.  The  second  electrode  consists  of  a 
lucite  plate  with  grooves  in  which  stainless  steel  rods  of  various  diameters 
can  be  placed.  These  rods  are  electrically  connected  to  the  sustainer 
capacitor.  An  optical  arrangement  of  two  flats  and  a 325  cm  radius  mirror 
(see  Fig.  10)  images  streamers  on  a photographic  plate  (Polaroid  High 
Speed  Type  57,  rating  ASA  3000),  arranged  so  as  to  produce  a slightly 
demagnified  image  with  magnification  m = 0.  92. 

Experiments  consist  of  filling  the  discharge  chamber  with  selected 
gas  and  firing  the  discharge,  photographing  the  streamers  generated  on  the 
rods  while  they  move  across  the  discharge. 

Since  we  do  not  have  time  discrimination  within  the  discharge 
pulse,  photographs  are  prepared  for  successively  longer  pulse  duration 
until  the  discharge  arcs.  When  this  happens  the  photograph  is  totally 
overexposed. 

1 . Eods  Parallel  to  the  Optic  Axis 

Three  rods  were  used  in  this  experiment,  placed  in  the  lucite 
electrode  holder  from  left  to  right  with  diameters  0.  16  cm,  0.  24  cm  and 
0.  32  cm  respectively,  and  using  it  as  cathode.  The  streamers  produced 
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for  discharges  of  different  duration  in  nitrogen,  using  applied  electric 
field  Eq  = 2600  V/cm  and  electron  density  n^^  = 10  cm  , are  shown 
in  Fig.  11;  the  rods  increase  in  size  from  left  to  right,  the  smallest  rod 
producing  the  largest  streamer  as  might  be  anticipated.  The  rectangular 
aperture  in  the  perforated  plate  is  placed  with  its  long  axis  normal  to  the 
optic  axis  so  as  to  locate  the  streamers  as  near  to  the  focal  plane  as 
possible,  while  the  rods  are  parallel  to  the  optic  axis.  It  is  immediately 
evident  that  the  streamers  from  the  thinnest  rod  have  propagated  most 
rapidly  across  the  discharge,  as  well  as  achieving  the  greatest  geometric 
complexity.  They  tend  to  follow  the  field  lines,  and  we  note  that  the  field 
distortion  at  the  rods  is  so  great  that  rod  position  is  more  important  than 
rod  radius  in  determining  streamer  growth  rate.  Nonetheless  we  can 
estimate  a growth  rate  from  Fig.  11;  we  get  a slowly  increasing  streamer 
velocity  of  ~ 1.5x10^  cm/sec. 

2.  Rods  Perpendicular  to  the  Optic  Axis 

Since  the  streamers  in  the  previous  case  are  following  rather 
complicated  field  lines,  their  velocity  is  not  very  well  defined,  and  a more 
uniform  field  would  be  desirable.  The  rectangular  EB  window  and  the 
cathode  rods  were  rotated  90°,  and  in  addition  the  two  outer  rods  were 
removed  and  the  lucite  electrode  holder  was  covered  with  aluminum  foil 
to  provide  a uniform  electric  field  in  the  discharge.  The  central  rod  pro- 
vided a discontinuity  on  which  streamers  were  likely  to  form,  and  the 
apparatus  was  adjusted  so  that  the  focal  plane  lay  along  the  central  rod. 

A typical  series  of  shots  is  shown  in  Fig.  12a.  These  were  per- 
formed in  nitrogen  at  atmospheric  pressure,  with  E^  = 2640  V/cm  and 
12  -3 

n^^  = 10  cm  . A halo  region  is  visible  around  the  streamer  tip,  as 
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Fig.  IZa  Streamer  Growth;  Single  Cathode  Rod  Perpendicular  to  Optic  Axis 


expected  in  the  model  discussed  above,  in  which  it  corresponds  to  the 
region  within  which  secondary  ionization  is  important.  The  streamer 
length  c is  plotted  vs  time  as  the  dotted  line  in  Fig.  12b,  in  which  it  can  be 
seen  that  the  streamer  is  accelerating  across  the  discharge  gap.  We  will 
discuss  this  further  below. 

The  intensity  of  the  image  enables  us  to  estimate  the  temperature 
of  the  streamer  plasma.  The  fluence  at  the  image  plane,  neglecting  losses, 
will  be: 


1. 

image 


2 2 
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TTbE 
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J/cm 
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where  a is  concave  mirror  stop  radius,  m is  magnification,  R is  distance 

from  source  to  concave  mirror,  E is  energy  radiated  from  streamer  in 
3 

joules  per  cm  , and  b is  a characteristic  dimension  of  the  streamer  (we 

take  the  radius).  We  note  that  the  fluence  required  to  produce  10%  increase 

118  2 

in  density  above  background  on  film  is  F = — — erg/cm  , where  S is 
the  ASA  number,  at  X = 5550  R.  A lower  limit  to  the  streamer  energy 
radiated,  which  completely  exposes  the  film,  is  then  - lOF,  which 

gives 


8m^  R^  1.  18  X 10"^  r / 

E = — ^ J/cc. 

7ra  b 


Taking  m = 0.94,  b = 0.  1 cm,  a = 10  cm,  R = 312  cm,  S = 3000,  we 

6 3 * 

get  E = 8.  6 X 10  j/cm  , or  for  exposure  ~ 5 p,sec,  at  tip,  E = 1. 72 

watt /cm 

(9) 

Using  the  tables  for  radiation  from  heated  air,  which  will  approxi- 
mate that  for  N^.  we  observe  that  this  radiation  corresponds  to  a streamer 


-41- 


VI 


temperature  T > 6000°K.  Taking  the  halo  as  ~ 10%  of  this  radiation  implies 
a halo  temperature  T ~ 4000°K.  We  note  that  the  Wien  blackbody  radiation 
peak  occurs  at  4800  R and  7244  R for  these  two  temperatures,  so  that  using 
the  above  estimate  for  fluence  on  the  film  is  approximately  self  consistent. 
Thus,  the  streamer  appears  to  have  a temperature  close  to  that  required  by 
the  consideration  of  heat  balance  in  the  foregoing  model. 

B.  FRAMING  CAMERA  EXPERIMENTS 

Recently,  we  have  obtained  much  more  precise  data  about  streamers 
in  our  discharge  by  using  a high  speed  framing  camera  to  record  successive 
stages  in  the  development  of  a single  streamer.  This  camera,  which  has  a 
maximum  framing  rate  of  44,  000  frames  per  second  is  shown  in  Fig.  13 
mounted  to  view  the  discharge  through  a window  in  the  bottom  of  the  cell. 

The  camera  was  aligned  to  photograph  a narrow  strip  of  discharge  extend- 
ing from  the  cathode  to  anode.  To  ensure  that  the  largest  streamer  would 
form  in  the  region  observed,  a loop  of  fine  tungsten  wire  was  affixed  to  the 
cathode  at  that  point.  In  this  way,  we  were  able  to  observe  the  development 
of  cathode  streamers,  as  can  be  seen  from  Fig.  14,  which  shows  a typical 
framing  sequence. 

Data  obtained  from  successive  photographs  of  a single  streamer  is 
of  course  more  accurate  than  that  obtained  from  integrated  photos  of  suc- 
cessive shots,  and  the  framing  camera  also  permits  us  to  photograph 
streamers  which  lead  to  arcing.  (In  Fig.  14,  the  last  few  frames  were 
slightly  fogged  by  the  subsequent  arc.  ) The  framing  camera  has  also 
greatly  increased  our  rate  of  data  acquisition,  since  one  shot  now  provides 
as  much  data  as  did  many  shots  under  the  old  procedure. 


i 
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HIGH  VOLTAGE  AREA 


FRAMING  CAMERA 


The  streamer  photographs  in  Fig.  14  agree  quite  well  with  the  model 
in  Fig.  4.  The  streamer  does  indeed  have  the  form  of  a cylinder  with  a 
ball  on  the  end.  (The  ball  is  barely  visible  in  Fig.  14;  it  is  more  clearly 
visible  in  other  framing  sequences.  The  ball  is  always  seen  when  the  ex- 
posure is  sufficient.  ) The  framing  photographs  were  the  first  experimental 
confirmation  of  these  predicted  details  of  the  streamer  structure,  and  the 
evident  agreement  with  the  model  is  a significant  confirmation  of  the  assump- 
tions. 

Figure  15  is  a plot  of  the  data  in  Fig.  14.  As  can  be  seen,  within 
the  accuracy  of  these  measurements,  the  streamer  growth  is  exponential. 

The  occurrence  of  an  arc  corresponds  to  the  streamer  reaching  a length  of 
about  1.  5 cm.  (The  anode  to  cathode  distance  is  one  inch.  ) This  behavior, 
which  is  quite  reproducible,  may  indicate  that  ionization  extends  beyond 
the  visible  streamer,  but  is  more  likely  due  to  the  onset  of  a faster  in- 
stability. The  streamer  is  a slow  mode  which  occurs  in  discharges  below 
the  conventional  arcing  threshold.  When  the  streamer  has  propagated  more 
than  halfway  across  the  medium,  a very  strong  field  is  produced  between 
the  streamer  tip  and  the  anode,  exceeding  the  threshold  for  rapid  break- 
down of  the  remaining  segment  of  the  path  needed  for  an  arc. 

The  streamer  appears  to  start  with  a length  of  about  1 millimeter. 

In  Fig.  14,  this  may  merely  represent  the  size  of  the  wire  loop  on  the 
cathode,  but  since  the  arcing  time  in  this  discharge  is  not  greatly  changed 
by  the  addition  of  the  wire,  streamers  of  this  size  apparently  form  rapidly 
at  the  cathode  in  any  case. 
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Fig.  15  Growth  of  the  Streamer  Shown  in  Fig.  14 
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Data  like  those  in  Fig.  14  wer'e  taken  over  a wide  range  of  discharge 
currents  and  electric  fields  (but,  only  in  nitrogen  at  atmospheric  pressure.  ) 
We  find  that  streamer  growth  is  always  nearly  exponential,  so  each  shot  can 
be  described  by  an  observed  growth  rate.  To  examine  the  dependence  upon 
discharge  conditions,  the  data  were  grouped  into  slow,  intermediate  and  fast 
growth  cases.  The  distribution  of  these  points  in  the  current-power  plane  is 
shown  in  Fig.  16. 

A clear  dependence  is  apparent.  Although  the  points  cover  a range 
of  more  than  threefold  variation  in  growth  rate  (which  represents  a very 
large  difference  in  streamer  size,  since  this  is  exponential  growth,  which 
we  typically  observed  over  several  E-folding  times),  they  all  lie  in  a fairly 
narrow  strip.  This  region  in  P - j is  the  threshold  of  the  streamer  insta- 
bility. Indeed,  the  straight  line  shown  in  Fig.  16  cleanly  separates  the 
slow  growth  cases  from  the  intermediate  and  fast  growth  cases.  This  line, 
which  is  quite  precisely  determined  by  the  data,  corresponds  to 

P - 3.  8 j = 1 

3 2 

where  P is  in  kW/cm  and  j is  in  A/cm  . 

In  the  P - j plane,  a straight  line  through  the  origin  would  represent 
a constant  value  of  electric  field.  The  observed  instability  line  is  almost, 
but  not  quite  of  this  form.  The  difference  is  apparent  when  the  data  are 
plotted  against  electric  field.  This  is  done  in  Fig.  17,  where  the  same  data 
points,  and  the  same  separation  line  are  shown  in  the  P - E plane.  Here 


The  electric  field  is  set  by  the  applied  voltage;  the  current  is  determined 
by  this  and  by  the  e-beam  current,  which  can  be  independently  selected. 
Thus,  there  are  two  electrical  parameters,  which  can  be  variously  stated, 
e.  g.  , current  and  electric  field,  or  current  and  power  density,  etc. 


P (kW/cm3) 


the  line  is  a hyperbola.  At  large  powers,  the  line  is  nearly  vertical,  im- 
plying that  y is  there  dependent  primarily  upon  E.  But  at  lower  powers 
and  higher  electric  fields,  the  line  is  nearly  horizontal,  indicating  that  y 
is  mainly  P dependent  in  that  region. 


C.  DISCUSSION  OF  RESULTS 
1.  Comparison  with  Theory 

Using  Eq.  (18)  we  can  calculate  streamer  length  versus  time  for  the 

nitrogen  discharge  shown  in  Figs.  12a  and  12b,  for  which  E^  = 3000  V/cm, 

“4  o _ 

cr^  — 2.  7x10  ,N  = 1 and  T = 6200  K.  The  power  loading  is  thus  approx- 
imately  2.  4 kW/cm  ; there  is  some  uncertainty  in  since  electron  density 
has  not  been  directly  measured,  due  to  the  electrode  geometry.  The  streamer 
lengths  calculated  for  various  power  loadings  are  shown  in  Fig,  12b,  together 
with  the  measured  experimental  curve.  For  most  of  its  existence  in  this  dis- 
charge the  streamer  should  be  moving  at  sonic  velocity.  There  appears  to  be 
farily  good  agreement  between  experiment  and  theory  in  this  case.  Consid- 
erable data  also  exists  on  arcing  times  in  discharges  of  various  dimensions. 
Since  we  can  approximately  identify  the  arcing  time  with  the  streamer  transit 
time  across  the  discharge,  a comparison  between  these  results  would  be  in- 
teresting. In  Fig.  18  the  expected  transit  time  for  streamers  of  various  lengths 
in  the  atmospheric  pressure  gas  mixture  He;  N^:  CO^  3:2:1  is  shown  as  func- 
tion of  power  loading.  Experimental  measurements  of  arcing  limits  for  the 
discharge  device  used  in  the  above  experiments,  with  electrode  spacing 
D = 2.  54  cm,  are  shown  as  crosses  and  open  squares.  One  experimental 
arc  was  obtained  in  another  device  with  spacing  D = 20  cm,  which  is  shown 
as  a solid  circle.  The  experimental  scatter  may  be  caused  by  streamers 
simultaneously  initiating  on  opposite  sides  and  meeting  halfway,  thus  reduc- 
ing the  arcing  time.  There  is  also  some  uncertainty  in  the  delay  time  before 
streamer  initiation,  which  would  have  the  opposite  effect. 

Experiments  have  also  been  performed  in  a device  with  D=4  cm  which 
can  be  operated  at  200°K  or  300°K.  The  effect  of  starting  at  low  temperature 
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will  be  greatest  during  the  sonic  propagation  regime,  since  the  streamer  will 
move  at  a lower  velocity  during  this  phase.  T’iC  subsonic  phase  will  change 
very  little,  however,  since  the  energy  required  to  go  from  200°K  to  300°K  is 
a very  small  fraction  of  that  required  to  reach  6200°K.  In  the  supersonic 
phase  the  streamer  will  move  at  2/3  its  previous  velocity  due  to  the  increased 
gas  density.  Experimental  results  for  arcing  times  together  with  theoretical 
streamer  transit  times  are  shown  in  Fig.  19.  The  gas  mixture  is  85%  N2,  10% 
CO_  and  5%  CO,  so  that  its  thermal  properties  are  dominated  by  N_.  The  pres- 
ence  of  carbon  atoms  in  the  heated  streamer  gas  implies  that  the  stable  streamer 
temperature  will  be  near  5350°K  (see  Section  II,  Fig.  7),  since  the  ionization 
potential  is  lower.  Again,  the  experimental  scatter  may  be  due  to  simultaneous 
formation  of  streamers  which  then  meet  in  the  middle.  Although  the  200°K 
cases  are  close  to  the  4 cm  discharge  theoretical  time,  the  300°K  cases  arc 
out  some  10%  sooner  than  predicted.  This  is  probably  caused  by  longitudinal 
electron  beam  shields  present  at  the  edges  of  the  discharge. 

Since  the  framing  camera  studies  were  done  only  for  pure  nitrogen 
discharges  at  one  atmosphere,  that  data  does  not  permit  a comparison  with 
the  predicted  pressure  dependence  or  ionization  potential  scaling.  However, 
we  can  test  the  theory  against  the  observed  form  of  the  streamers,  their 
temporal  growth,  and  the  scaling  of  growth  rate  with  discharge  parameters  dis- 
played in  Figs.  16  and  17.  As  already  noted,  the  framing  camera  photographs 
do  show  a streamer  structure  composed  of  a luminous  column  of  gas  with  a 
ball  on  the  end,  as  sketched  in  Fig.  4. 
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treamer  Transit  Time  for  Cavities  with  Various  Separations 
redicted  for  200°K  and  300°K  N2  Discharge  at  1 atm.  Experi- 
lental  arcing  times  are  shown  for  200°K  and  300°K  in  N7:  CO7:  CO 
.84:0.01:0,06. 
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The  framing  camera  data  also  show  the  streamers  growing  exponen- 
tially. Such  growth  is  quite  consistent  with  the  theory,  as  can  be  seen  by 
rewriting  Eq.  (15)  for  the  subsonic  velocity  as 

2 

^^0^0  rai  r? 

dt=  A.keN  [^[4 

If  c/a,  the  streamer  aspect  ratio,  is  constant  during  growth  , this  can  be 
written 

where  the  growth  rate  is 

3a  E 

o 

^ ^ X k e 

Hence,  unless  c/a  changes  strongly  with  discharge  parameters,  the  theory 

2 

predicts  a growth  rate  proportional  to  the  power  loading,  E^  . From 
Fig,  17,  one  sees  that  this  is  true  at  high  electric  field,  but  not  in  the  weak- 
field,  high-current  regime.  Since  the  assumptions  underlying  our  theory 
(e.g.  , that  the  corona  around  the  streamer  tip  extends  well  beyond  the  radius 
of  the  streamer  column)  are  appropriate  to  the  strong  field  operating  mode, 
one  expects  agreement  with  experiment  there,  as  is  borne  out  by  our  data. 

In  the  high  current,  low  field  regime,  the  resistive  heating  of  the  gas  beyond 
the  streamer  tip,  i.e,,  beyond  the  region  of  the  breakdown  field,  would  have 
to  be  considered.  Since  such  heating  has  not  been  considered  in  the  present 
analysis,  the  lack  of  agreement  there  is  not  surprising. 

This  conclusion  is  also  supported  by  the  lack  of  a sonic  barrier  to 
streamer  propagation  under  such  conditions.  Several  of  the  framing  camera 

From  the  theory,  one  expects  this  to  be  at  least  approximately  true, 
an  expectation  which  is  borne  out  by  the  framing  camera  photographs. 
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sequences  showed  exponential  growth  continuing  unimpeded  into  a velocity 
range  which  would  be  supersonic  for  the  background  gas.  Evidently,  the 
streamer  was  growing  into  gas  already  heated  by  the  intensified  current 
ahead  of  the  streamer.  Such  heating,  which  turns  on  smoothly  as  the  streamer 
approaches,  raises  the  sound  speed  and  hence  prevents  any  abrupt  change  in 
behavior  as  the  growth  goes  supersonic  relative  to  the  cold  background  gas. 
Again,  such  behavior  is  characteristic  of  the  current-dominated  regime, 
which  is  not  of  much  interest  for  a laser  amplifying  medium.  In  the  regime 
which  is  of  interest,  the  theory  agrees  well  with  our  data. 

2.  Predictions  of  Arcing  Limits  for  Discharges 

Streamer  behavior  is  most  concisely  described  by  calculating  a maxi- 
mum loading  energy  E = i E t..  where  r.  is  the  time  taken  for  stream- 
ers  to  cross  the  discharge.  Using  the  appropriate  values  from  Eq.  (18),  these 
loading  energies  have  been  calculated  for  the  gas  mixture  H e:  N^:  CO^:  3:  Z:  1 at 
atmospheric  pressure,  and  are  shown  in  Fig.  20a  as  a function  of  the  dis- 
charge power  loading.  The  parameter  D is  the  discharge  electrode  spacing 
in  centimeters.  Physically  the  flat  portion  of  the  curve  at  low  power  corres- 
ponds to  the  subsonic  regime.  For  wider  discharges  the  sonic  and  super- 
sonic regimes  are  reached  in  the  discharge;  since  the  streamer  is  as  it  were 
delayed  by  the  sonic  regime  more  energy  can  be  introduced  before  the  streamer 
transits  and  causes  arcing.  Increasing  power  input  causes  the  sonic  regime 
to  move  closer  to  the  initiating  electrode.  An  experimental  point  correspond- 
ing to  the  only  observed  arc  in  a device  with  D = 20  cm  is  indicated. 

As  the  pressure  is  reduced  tlie  thermal  diffusivity  of  the  gas  increases. 
Consequently  the  streamer  radius  b,  and  the  diameter  of  the  cap  sphere  2a, 
will  both  increase  from  Eq.  (16).  This  results  in  relatively  slower  longitudinal 


movement  of  the  streamer,  since  while  the  current  collected  increases  as 
the  square  of  the  radius,  the  gas  to  be  heated  increases  as  the  cube,  causing 
slower  heating.  Consequently,  at  low  pressures  more  energy  can  be  intro- 
duced per  molecule  of  gas  than  at  high  pressures  before  the  streamer  reaches 
the  opposing  electrode  and  arcing  ensues.  The  predictions  for  the  same  gas 
at  l/lO  atm  are  shown  in  Fig.  20b.  We  note  that  it  has  been  known  empiri- 
cally for  some  time  that  arcing  occurs  "later”  in  EB  discharges  at  low  pres- 
sure than  at  high  pressure.  In  the  present  model  this  is  due  only  to  the  de- 
pendence of  the  streamer  lateral  dimensions  on  the  gas  thermal  diffusivity. 

It  is  interesting  to  see  that  for  typical  CW  discharge  conditions  we 
would  expect  to  be  almost  wholly  in  the  subsonic  regime.  Consequently  the 
discharge  energy  loading  will  be  lower  than  at  the  high  power  loading  con- 
ditions typical  of  pulsed  operation,  despite  the  effect  of  increased  diffusivity 
mentioned  above.  This  is  the  reason  for  the  lower  specific  energy  input 
possible  in  CW  discharges.  A point  corresponding  to  HPL-lO  (which  actually 
operates  in  He:  N2,:  CO2:  CO  6:  4;  1:  1 at  l/lO  atm)  is  indicated  in  Fig.  20b, 
which  was  calculated  for  3 : 2 : 1 : 0.  Since  these  two  gas  mixtures  have  very 
similar  thermal  properties  and  effective  ionization  potentials,  their  behaviour 
should  be  close.  In  fact  it  appears  that  HPL-10  arcs  at  a lower  loading  than 
predicted  on  this  model.  This  could  be  due  to  the  presence  of  low  ionization 
potential  contaminants  in  the  gas. 

3.  Effect  of  Impurities 

The  effect  of  a small  amount  of  easily  ionized  substance  on  the  dis- 
charge streamer  temperature  has  been  discussed  above.  We  observe  from 
Fig.  7 that  0.  1%  molar  of  sodium  decreases  the  streamer  steady  state  tem- 
perature to  T = 3500°K,  and  consequently  the  streamer  cam  move  more 
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rapidly  since  it  heats  the  gas  to  a lower  temperature:  numerically  the  arcing 
time  is  reduced  to  70%  of  its  value  in  pure  gas.  It  is  clear  that  any  potas- 
ium  compounds  present  will  have  an  even  more  drastic  effect.  The  presence 
of  easily- ionized  material  will  tlius  increase  the  likelihood  of  arcs,  and  for 
discharges  operated  just  below  the  arcing  limit,  as  is  usually  the  case,  we 
expect  even  very  small  amounts  of  contaminant  to  initiate  arcing.  In  CW  dis- 
charges the  same  effect  would  be  produced  by  any  low  ionization  potential 
organic  compound  present  in  the  gas  flow. 


D,  INITIATION  OF  STREAMERS 


The  process  of  streamer  initiation  at  electrodes  usually  occurs 
on  a discontinuity  or  projection.  This  presumably  takes  place  because 
the  electric  field  intensification  near  a projection  results  in  more  rapid 
heating  of  the  gas  at  that  point,  and  this  continues  until  the  gas  reaches 
a temperature  at  which  its  conductivity  becomes  high  enough  to  begin 
propagating  as  a streamer,  as  discussed  previously. 

There  must  be  a critical  size  for  such  a projection  to  cause  a 
streamer,  for  if  the  projection  is  very  small  conductive  gas  cooling  will 
be  so  rapid  that  the  steady- state  temperature  reached  will  be  below  that 
required  for  a stable  streamer.  Although  streamers  could  ultimately 
form  at  such  projections,  this  would  be  on  a time  scale  corresponding 
to  that  required  to  heat  the  entire  discharge  to  a conductive  temperature 
rather  than  a small  part  of  it. 

We  will  estimate  the  heating  that  occurs  near  a spheriodal 
protrusion  on  an  electrode. 

Considering  a spheroidal  protrusion  of  height  c and  base  radius 

b,  we  can  write  the  field  enhancement  near  the  tip  as 

ft  -£ti£_ 


\L^  i\ . 
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]^o  - ij 
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where  qo  = " b^/c^)“^/^,  and  limiting  ourselves  to  a volume  whose 

radial  dimension  is  the  radius  of  curvature  of  the  spheroidal  tip,  which 

.2 

is  p (55  -5 — , we  equate  the  volumetric  conductive  cooling  from  Eq.  (1) 
with  the  volumetric  heating  rate,  and  obtain; 

A t7  t.2  „2 

^ P 


(20) 
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or 


2 A T ^ 

W y..  pZe2 

We  note  that  for  one- dimensional  planar  wall- dominated  conductive 
cooling  A a 2;  for  cooling  of  a cylinder  A « 15  as  before,  and  we  will 
use  the  latter  value. 

The  critical  base  radius  for  a spheroidal  protrusion  in  a nitrogen 

discharge  at  1 atm  is  shown  in  Fig.  21  as  fxinction  of  spheroidal  aspect 

2 

ratio  b/c  and  discharge  power  loading  a E^.  The  final  temperature  is 
taken  as  T = 6200°K.  For  b/c  = 1 we  have  a hemispherical  boss  on 
the  electrode,  and  for  this  the  maximum  field  is  always  three  times 
the  discharge  field,  regardless  of  radius,  so  that  the  heating  rate  at 
its  tip  is  always  nine  times  the  ambient  heating  rate.  However,  if  the 
initial  heating  takes  place  in  a thin  layer,  equating  layer  thicknesses 
with  radius  of  curvature  as  above  is  incorrect.  We  would  then  expect 
the  initiation  time  lag  to  be  much  greater  than  estimated  here,  and  to 
correspond  to  the  time  required  to  heat  up  the  electrode  surface  itself  — 
essentially  infinite  by  comparison  to  the  times  considered  here. 

If  the  protrusion  is  larger  than  the  critical  diameter  2b,  then 
volume  heating  will  take  place  and  the  effect  of  surface  heat  loss  will 
reduce  the  volume  heating  rate  only  slightly,  and  may  be  neglected. 

We  may  then  estimate  the  latent  time  taken  to  heat  the  gas  in  the 
neighborhood  of  the  protrusion  as: 


^ g (T) 

- 

where  e (T)  is  the  enthalpy  (see  Appendix  A)  required  to  raise 

unit  volume  of  gas  to  the  temperature  T at  which  the  electrical 
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DISCHARGE  POWER  LOADING  w/cm 


Fig.  21  Streamer  Critical  Radius  and  Latent  Period  for  Spheroidal  Protrusion 


conductivity  is  high  enough  for  streamer  propagation,  and  [j  E]  is  the 

power  loading  per  unit  volume  of  gas  near  the  electrode. 

The  electrode  power  loading  per  unit  volume  is  very  much  greater 

than  in  the  main  discharge,  and  especially  at  the  cathode.  As  pointed  out 

by  Cobine,  the  current  density  required  near  the  cathode  is  typically 
_2 

near  15  amp  cm  in  the  glow  discharge  at  high  pressure:  the  discharge 
contracts  so  that  electrons  actually  leave  the  cathode  only  over  a small 
fraction  of  its  surface.  For  a 500  V cathode  drop,  with  cathode  layer 

_3 

thickness  ~ 10  cm,  the  power  loading  will  be  3000  MW/cc.  Under  these 
conditions  T is  negligible  and  we  may  assume  that  streamers  form 
instantly  at  the  cathode  and  begin  their  progress  across  the  discharge. 

At  the  anode  there  will  be  some  delay  and  we  anticipate  preferential 
streamer  formation  on  discontinuities;  we  also  expect  them  to  be  shorter. 
Since  arcing  probably  corresponds  to  cathode  streamers  completing  the 
discharge  transit,  we  would  expect  the  transit  durations  discussed  previ- 
ously to  be  the  arcing  times,  and  would  predict  negligible  initiation  delay. 
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E.  CONCLUSION 


Simple  physical  assumptions  lead  to  an  analytic  model  for  streamers 
in  non- self- sustaining  discharges.  This  model  predicts  streamer  behavior 
similar  to  that  observed,  and  discharge  arcing  times  calculated  from  it  are 
quite  close  to  experimental  values,  for  a variety  of  pulsed  and  CW  discharges 
Framing  camera  photographs  have  confirmed  the  anticipated  structure  of  the 
streamer.  Scaling  laws  derived  for  streamers  are  in  accordance  with  ob- 
served results.  The  effect  of  low  ionization  potential  contaminants  is  to  re- 
duce streamer  temperature  and  thus  result  in  more  rapid  streamer  transit 
and  arcing.  This  model  will  be  developed  in  more  detail  in  the  future. 


APPENDIX  A 


GAS  SPECIFIC  HEATS 


A.  CONSTANT  VOLUME  CONSTANT  PRESSURE  (CVCP)  ENTHALPY 

In  the  subsonic  regime,  the  energy  required  to  raise  one  gram  of 
gas  at  constant  pressure  from  temperature  T^  to  T is  equal  to  the 
enthalpy  difference  between  these  temperatures.  Writing  enthalpy  as 
q(T),  the  energy  is  AQ  = q{T)  - q(T^)  joules  per  gram.  The  energy 
required  to  raise  unit  volume  of  gas  over  the  same  temperature  range 
at  constant  pressure  will  then  be 


= /.(!!) 


dT  . 


where  p is  gas  density  at  temperature  t' , and  C^  = . 


Under  these 


conditions  p = P^T^/t',  where  p^  is  density  at  T = T^,  and 


= Po^o  /'(It)  • 


Fortunately  e can  generally  be  represented  fairly  well  by  a simple 

power  law  for  most  molecular  gases.  We  note  that  in  the  case  of  a 

perfect  monatomic  gas  we  can  immediately  write 

€ = 5/2  N kT  fn  (T/T  ), 
o o ' o 

where  N is  initial  number  density  at  T = T and  k is  Boltzmann' s 
o o 

constant.  In  the  case  of  molecular  gases  the  energy  used  in  exciting 
different  vibrational  modes,  rotation,  and  dissociation  must  be  included, 


I 


as  well  as  the  fact  that  the  gas  density  under  constant  pressure  will  de- 
crease more  rapidly  than  1/T. 

1,  Helium 

Treating  this  as  a perfect  monatomic  gas  and  ignoring  ionization 

for  T < 15,00C"^K,  we  take  initial  conditions  as  T = 300°K  and  P = 1 atm: 
’ ’ o 


= 0.  254  fn  (T/T^)  j/cm  . 


(A-2) 


2.  Nitrogen 


Using  the  thermodynamic  tables  for  N21  Eq,  (A-1)  is  plotted  in 


Fig.  A-1.  A suitable  fit  for  the  region  of  interest  is 

€n  = €,(T/Tj)“, 

with  c = 1 j/cc  and  Tj  = 3500°K,  a = 0.  63. 

3.  Carbon  Dioxide 


(A-3) 


(11) 


The  constant  pressure  enthalpy  for  CO2  is  shown'  ' in  Fig.  A-2 

3 

The  best  fit  for  the  energy  required  to  heat  1 cm  to  temperature  T is  then: 


^C02  = J 


(A-4) 


where  T^  = 691°K  and  fj,  = 0.  71. 

B.  NUMERICAL  STREAMER  ENERGIES  IN  THE  SUBSONIC  REGIME 

Using  the  above  values  of  the  constant -volume  constant-pressure 
(C^Cp)  enthalpy  for  the  gas  mixtures  He:N2:C02'.C0  0:1:0;0,  3:2: 1:0, 
0:3:1:0,  16:3:1:0,  we  obtain  the  following  limiting  energies  at  which 
streamers  can  exist  in  these  gases  at  atmospheric  pressure. 


1.  Nitrogen 


€(T)  = (T/3500)^' N = 1,  T = 6200°K. 


• 68- 


Hence 


€ = 1.43  j/cm^ 

2.  He:N.,:CO.,:3:2:l 

The  most  readily  ionized  component  of  this  mixture  when  dissoci- 
ated will  be  the  10%  carbon  atoms,  and  using  this  in  the  Saha  equation  and 
solving  (5)  gives  the  line  shown  in  Fig.  7,  from  which  we  obtain  T ~ 5350°K 
for  streamer  formation  under  typical  conditions.  Using  the  CVCP  enthalpy 
for  each  component,  we  get 

€ w 1/6  [3  X 0.  254  fn  (T/300)  + 2 (T/3500)®- + (T/691)°’ , 
or 

€ « 1.  51  j/cc  at  5350°K. 

3.  He;N^;COo;CO  3;2;i;i 

This  gas  may  be  expected  to  have  virtually  the  same  thermal  be- 
havior as  3:2: 1:0  at  1 atm. 

4.  N^:CO^:Ht  3:1:0.  08 

We  may  apply  the  above  reasoning  and  obtain  a molar  fraction  11% 
carbon  atoms  in  the  dissociated  gas,  giving  essentially  the  same  critical 
temperature  T - 5350°K.  The  CVCP  enthalpy  is  then 

e « A [3  X {T/3500)°*  + (T/691)°' . 

or 

€ « 2.  05  j/cc  at  5350°K. 

5.  He:N2;C02t  16:8:1 

The  critical  temperature  in  this  gas  will  be  near  6000°K  due  to 
its  lower  carbon  content.  Using  Appendix  A gives  the  corresponding 
CVCP  enthalpy  as  6 »>  1.  12  j/cc  at  6000°K. 


C.  INTERNAL  ENERGY 


In  the  supersonic  regime  expansion  of  the  gas  during  the  heating 
process  may  be  neglected,  and  the  total  energy  required  is  the  difference 
in  internal  energy  between  the  lower  and  higher  temperatures.  Again,  for 
a monatomic  gas  with  negligible  ionization  this  will  be: 

e = I R A T 

while  for  molecular  gases  the  contributions  from  dissociation  and  ionization 
must  be  included. 

1.  Helium 

e = kATJ  cm’^ 

d o 

2.  Nitrogen 

In  the  thermal  regime  5000-6000°K,  for  gas,  initially  at  atmospheric 
pressure  and  300°K,  we  have: 

Thus  at  6200°K,  e = 13.  1 J cm'^  and  5350°K,  € = 11.3. 

3.  Carbon  Dioxide 

In  the  thermal  regime  5000-6000°K,  for  gas  initially  at  atmospheric 
pressure  and  300°K,  we  have: 

Thus  at  5200°K,  € = 13.  68  J cm"^  and  at  5350°K,  e = 11.  8. 

4.  He;N^;CO^  3:2:1 

For  this  mixture  initially  at  atmospheric  pressure  and  300OK,  the 
weighted  average  of  the  above  quantities  yields: 

€ = 1.  33  X 10"^  T J cm"^,  where  5000  < T s 6000. 

Thus  for  T > 53500K,  e = 7.  1. 
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